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Abstract. We review some of the recent progress in our understanding of the physics 
of ultrarelativistic heavy ion collisions due to applications of AdS/CFT correspondence. 
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1. Introduction 

^ I Ultrarelativistic heavy ion collision experiments are being carried out with the goal of 

>■ ■ creating a thermalized medium made out of deconfined quarks and gluons, known as the 

o 

I those characterizing a static thermal bath, such as the equation of state and the order 

^ I of the phase transition, to the dynamical properties, like the transport coefficients, 

including shear and bulk viscosities of QCD matter. With heavy ion experiments 
running at RHIC and LHC, qualitative and quantitative theoretical progress is required 
both to better understand the existing data and to suggest new observables to measure. 
_ Heavy ion collision is a multi-scale process, with the harder momenta dominating 

\ the early times right after the collision, and the softer momentum scales becoming 

^ ■ relevant at the later time scales. Due to the running of the QCD coupling constant 

the hard modes are weakly-coupled, while the soft modes are coupled strongly. In 
recent years a body of evidence has accumulated suggesting that the non-perturbative 
large-coupling effects could be more important for high-energy heavy ion collisions than 
was previously thought. Perhaps the most important piece of evidence is the success 
of ideal hydrodynamics in describing the evolution of the medium produced in heavy 
ion collisions [H [5] . Corrections to ideal hydrodynamics are inversely proportional to 
the coupling constant and are small when the coupling is large. An example of this is 
the shear viscosity, which (at weak coupling) scales as ~ / ln{l / g) in a thermal 
QCD medium with the temperature T |6]: extrapolating this result one sees that the 
shear viscosity is small when the coupling g is large, generating a small correction to 
the energy-momentum tensor of ideal hydrodynamics. Another argument in favor of 
strong-coupling djTiamics playing an important role in heavy ion collisions is the large 
amount of jet quenching needed to describe the data coming from RHIC and LHC. The 
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relevant jet quenching parameter g, which is proportional to the energy loss per unit 
path length of an energetic parton traversing a medium [7] , was found to be rather large 
in some analyses, possibly indicating that non-perturbative strongly-coupled effects are 
at work [8]. Note however that q has mass dimension of M^, and the existence of 
a large momentum scale q^^^ likely indicates that at least a part of the process is in 
fact weakly coupled, and, hence, perturbative. The last piece of evidence suggesting 
possible importance of strong coupling effects in the QCD plasma comes from lattice 
simulations, which indicate that the energy density and pressure of the thermal QCD 
medium is about 80 -^ 85% of that in an ideal gas of quarks and gluons for a broad 
range of temperatures above the temperature of the QCD phase transition Tc [HI ITO] . 
As the strong coupling constant as{T) should get small at large-T, one expects the 
high-temperature QCD medium to approach the Stefan-Boltzmann ideal gas behavior 
as T increases: the approach does take place, but is much slower than expected by in the 
perturbative QCD calculations, possibly indicating that non-perturbative effects are still 
important even at reasonably high-T. Note however that the resummed perturbation 
series may successfully describe the lattice data pT| . 

While none of the arguments presented above definitively demonstrates dominance 
of strong-coupling effects in heavy ion collisions at RHIC and LHC, we will proceed 
under the assumption that the strong-coupling effects are important and try to analyze 
heavy ion collisions assuming that the coupling is large. One reason for this is that 
exploring heavy ion collisions in the strong-coupling regime is very interesting by itself, 
providing a new theoretical angle on the process. In addition, it may be that the typical 
strong coupling constant in these collisions is neither very small nor very large, such 
that both small- and large-coupling approaches would have some chance of describing 
the data. 

The strong-coupling QCD analysis of such an ultrarelativistic process as a heavy ion 
collision is impossible with the present state of the QCD theory. Instead an opportunity 
to study the strongly-coupled field theories came up due to advances in string theory, 
culminating with the formulation of the Anti-de Sitter space/ Conformal Field Theory 
(AdS/CFT) correspondence [12[ [131 [HI [151 [16] , which is the duality between the J\f = 4 
SU{Nc) super Yang-Mills (SYM) theory in four flat space-time dimensions (conformal 
field theory, or simply CFT) and the type IIB superstring theory on AdSsxS^. In 
the limit of large number of colors A^^^ and large 't Hooft coupling A = g'^N^ such 
that A'c ^ A ^ 1, AdS/CFT correspondence reduces to the gauge-gravity duality: 
A/" = 4 SU(Ac) SYM theory at A'c > A > 1 is dual to the (weakly coupled) classical 
supergravity in AdSsxS^. The AdS/CFT correspondence gives us a powerful tool one 
can use to systematically study a strongly-coupled gauge theory (specifically. A/" = 4 
SYM theory). Using the AdS/CFT correspondence one can try to study heavy ion 
collisions at strong coupling. Indeed one should be mindful about the many differences 
between QCD and the N' = 4 SYM theory: the former is confining, while the latter is 
not; QCD coupling runs, while it is a constant in the conformal M = 4 SYM theory; 
QCD plasma has a chiral phase transition, while Af = 4 SYM plasma does not. Until the 
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errors introduced by using a theory different from QCD are understood, any conclusions 
we derive from applying the AdS/CFT correspondence to the study of QCD-mediated 
processes should be taken as qualitative at best. However, in some cases one may 
hope that AdS/CFT predictions are universally valid for all strongly-coupled theories, 
including QCD. This hypothesis is supported by some numerical successes of AdS/CFT, 
an example of one of which will be shown below. 

2. Elements of AdS/CFT correspondence and the strongly-coupled 
supersymmetric plasma 

For reviews of the AdS / CFT correspondence we refer the reader to [lEl [TTl [181 HH] • Since 
our goals in this mini-review are rather applied, we will simply state the tools we need 
to accomplish them, without presenting any proofs, and illustrate these statements by 
analyzing some of the properties of the supersymmetric plasma. 

AdSs X S^ space is obtained by stacking parallel D3 branes on top of each other. 
The resulting metric for (the Poincare wedge of) AdSsxS^ in the z ^ L limit is 
L2 



-dr + dx'^ + dz" 



+ dQi (1) 



where (t, x) = (t the regular four-dimensional coordinates, z G [0, +oo) 

is the fifth dimension of AdSs, and the branes are located at z = oo. L is the radius of 
the S^ sphere. For our purposes it is important to note that the AdSs metric in Eq. ([1]) 
(everything on the right except for the last term) satisfies Einstein equations 

Rmn — - Qmn -R + A qmn = (2) 

with M, = 0, . . . , 4 and the cosmological constant A = — 6/L^. 

The AdS/CFT correspondence is the statement of equivalence of the J\f = 4 SU(A'c) 
SYM theory in four space-time dimensions and the type IIB superstring theory on 
AdSsxS^ [l2]. For practical calculations the duality means the following. Each operator 
O in the four-dimensional CFT is paired up with the dual field if in the AdSs bulk. Then 
the duality can be mathematically formulated as [131 El 115] : 



exp 



d'^x 0{x) ^po{x) 



CFT, ^^^"^g 



(p{x,z] 



z'^ ^^q{x) 



.(3) 



On the left hand side of Eq. (I3|) we have an expectation value taken in the full four- 
dimensional CFT, while on the right there is a partition function for the type IIB 
superstring theory, calculated with a particular boundary condition on the bulk field 
^{x^z). A is the scaling dimension of the operator O. One can show that in the 
A'c ^ A ^ 1 limit the string partition function becomes classical, ^gtring ~ e^^^M 
with Si,uik the classical action on AdSsxS^ calculated with the ^{x^ A\zk.q ~ V^o(2;) 
boundary condition. Thus the incredibly complicated quantum field theoretical problem 
of calculating an expectation value in the full CFT is reduced to a classical calculation 
on AdSsxS^. This is a great simplification due to the AdS/CFT correspondence. 
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As an important example of the use of AdS/CFT duality which is also useful for 
heavy ion collisions, let us consider a thermal medium made out of strongly-coupled 
A/" = 4 SYM matter. Due to thermodynamic properties of black holes, the dual of the 
thermal medium in the gauge theory is the AdS Schwarzschild black hole (AdSSBH) 
with the metric [12] 



1 - 



dt^ + dx^ + 



dz^ 



z' 



h, 



1 - 



(4) 



The horizon of AdSSBH metric is located ai z = Zh- Introducing the Euclidean time 
tE = it the metric (jl]) can be extended to the Euclidean space-time. Expand the 



Euclidean version of the metric (jl]) near the horizon by redefining z 
the lowest non-trivial order in p we get 
4p2 



ds' 



dt% 



l2 

— ^ dx"^ + dp^. 

4 



at 



(5) 



We see that the (t^;, p)-part of the metric looks like the polar coordinates in 2d, with the 
metric ds2d = p^ dcp"^ + dp'^, if we define the "angle" (p = 2tE/zh- The metric ds2d has no 
kink singularity at p = only if it is periodic in (j) with the period of 2 vr. This means that 
the Euclidean time tE should be periodic with the period of AtE = n Zh- Remembering 
that a thermal field theory in Euclidean space is periodic with the period At^; = /3 = 1/T 
with T the temperature, we (i) find the Hawking temperature Th = l/ij^Zh) of the AdS 
Schwarzschild black hole (jl]) and (ii) identifying the Hawking temperature of AdSSBH 
with the temperature of the dual gauge theory we obtain 

ZH = (6) 

The metric of (jl]) is illustrated in Fig. [H We see from Eq. (j6]) that at high temperature 
T the horizon is close to the boundary of the AdS space at z = 0, while for low T the 
horizon is deep in the AdS5 bulk. 








(t,x) 



Figure 1. Pictorial illustration of the AdS Schwarzschild black hole metric (|4]). 
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Identifying the AdS Schwarzschild black hole as the gravity dual of the strongly- 
couples SYM plasma allows us to calculate the entropy of the SYM plasma by equating 
it with the Bekenstein-Hawking entropy of a black hole 

Sbh = (7) 

where A is the horizon area and Gn is the Newton's constant. Since the area of a 
5-sphere is A^s = L^, the horizon area for AdSSBH can be calculated via 

A = j dx^ dx^ dx^ ^Jdetg^^\^^ A55 (8) 

where det g^d = L^/z^ is the determinant of the part of the metric (jlj) in three dimensions 
{x^,x'^,x^). Using Eq. (EJ we get A = Vs^tt^T^ with V^d the 3d volume occupied by 
the plasma. Using this area in Eq. ([7]) along with 

O. - (9) 

we obtain the Bekenstein-Hawking entropy of AdSSBH, or, by AdS/CFT duality, the 
entropy of the strongly-coupled Af = 4 SYM plasma [20] 



S=^N^^T'V,d. (10) 

We have derived an expression for the entropy of a strongly-coupled medium in a thermal 
field theory by doing a very simple gravity calculation! Indeed the strong-coupling 
expression ffTOj) corresponds to a sum of an infinite number of Feynman diagrams: 
AdS/CFT duality allowed us to arrive at the answer in a much simpler way. 

Eq. (ITOi) should be compared to the Stefan-Boltzmann entropy of an ideal gas of 
non-interacting particles in A/" = 4 SYM theory 

SsB = ^N!T'Vsd. (11) 

We see that S/Ssb = 3/4: the entropy of a strongly-coupled supersymmetric plasma 
is 3/4 of that for the ideal gas. As was mentioned before, this result seems to be 
in a semi-quantitative agreement with the lattice QCD calculations for entropy above 
Tc giving S/SsB = 80 85% [9l [10], as illustrated in Fig. [2] for the entropy density 
s = S/V^d- Note, again, that in interpreting this semi-agreement one has to remember 
that M = A SYM and QCD are different theories, and the AT = 4 SYM prediction ([TOD 
for S/T^ would correspond to a horizontal straight line in the plot of Fig. [2] exhibiting 
no temperature dependence. 

Now suppose we want to find the expectation value of the energy-momentum tensor 
(EMT) {T^y) in the A/" = 4 SYM plasma (/i, = 0, . . . , 3 label the 4d space-time). Since 
gravity couples to the EMT of the matter fields in 4d, by AdS/CFT rules ((Sj) the EMT 
operator is dual to the metric tensor qmn in the bulk [21]. The AdS/CFT duality ([3]) 
leads to the following prescription for the calculation of (T^^), known as the holographic 
renormalization [21] . First cast the metric in the Fefferman-Graham form [22] 

ds'^ = — [g^ydx^dx" + dz^\ . (12) 
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Figure 2. Entropy density s (scaled down by T^) of tlie (2+l)-flavor QCD as a 
function of temperature determined by the lattice simulations in [TU]. The ideal gas 
Stefan-Boltzmann limit is shown by a horizontal line. (Reprinted figure with permission 
from [TU]. Copyright 2009 by the American Physical Society.) 




Then expand the 4d metric g^y{x,z) near the boundary of AdSs at 2; = 0. One can 
show that if the expansion starts with the Minkowski metric rj^i, at 2 = (as is the case 
for AdSs) then the next term in the expansion is order-z^, such that [21] 

g^y{x, z) = r]f,y + g^^Jix) + .... (13) 

With the expansion ( IT3l) in hand one can find the expectation value of the EMT by the 
following simple relation 



(14) 



2 7r2 

Applying this prescription to the AdSSBH metric (j4]) we first perform a substitution 

z 



z 



(15) 



1 + F 



with Zh = Zh a/2 recasting Eq. (j4]) into Feffer man- Graham form [23] 

2 



l2 

z^ 



1 - 



1 + 



z* 



— \ dx + dz 



z2 



(16) 



Reading off the z'^ coefficients of g^v{x, z) from Eq. (fT6ll we obtain the EMT of a static 
supersymmetric plasma 



vr 



(rM<^) = ^^c^ diag{3,l,l,l}, 



(17) 



As we will see below, holographic renormalization is a powerful tool that can be used 
to study heavy ion collisions at strong coupling. 
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3. Heavy ion collisions in AdSs 

As we mentioned above, heavy ion collision is a multi-scale process involving an 
interplay of small- and large-coupling phenomena. Unfortunately no single theoretical 
framework incorporates both perturbative and non-perturbative effects. Therefore, 
applying AdS/CFT duality to heavy ion collisions one has to make a choice of either 
treating the whole collision as a strongly-coupled process, or, perhaps more realistically, 
limit the application of AdS/CFT duality to the description of the produced medium 
at the later times after the collision when the QCD coupling is more likely to be large. 
In this Section we will discuss the progress in the former direction, while the advances 
in the latter one are presented in the next Section. 

3.1. Shock wave collisions 

We begin by modeling the whole heavy ion collision in the strongly-coupled framework 
of the AdS/CFT correspondence. The goal here is rather ambitious, and includes 
understanding particle production in the collision, followed by the thermalization of the 
produced matter, along with the subsequent evolution of the resulting thermal medium: 
essentially one wants to understand the whole collision in a single unified framework. 
While there is no nuclei (or other bound states) that we could imagine colliding in the 
conformal A/" = 4 SYM theory, one can follow Bjorken [25j and model the colliding 
nuclei as the Lorentz-contracted thin sheets of matter, which are uniform and extend 
to infinity in the transverse direction. The collision of two such "nuclei" is illustrated 
in Fig. [3] using light-cone coordinates = (x 




nucleus 1 nucleus 2 



Figure 3. The space-time picture of the ultrarelativistic heavy ion collision in the 
center-of-mass frame. The collision axis is labeled x^, the time is a:°. 

The nucleus moving along the x"'"-axis ("nucleus 1" in Fig. [3]) has a large T^"*" 
component of its EMT, with all other EMT components negligibly small. Moreover, 
due to assumptions we made about the uniform distribution in the transverse plane, 

is a function of x~ only with its support localized near x~ = due to Lorentz 
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contraction. The metric dual to such an ultrarelativistic nucleus which solved Einstein 
equations ([2]) is that of a gravity shock wave in the bulk |23] 

ds'' = ^ l-2dx+dx- + ^{T++{x~))zUx-^ + dxl + dzA (18) 
[ J 

with dx^ = (dx^)"^ + {dx"^)"^. The metric dual to "nucleus 2" in Fig. [3] is obtained from 
Eq. (ITS]) by the + ^ — interchange. 

The problem of understanding heavy ion collisions in this version of AdS/CFT 
approach can be formulated as follows: we know the metric for x"*" < or for x^ < 0, 
where it is simply a sum of the two shock wave metrics [Eq. (fTSl) and the one obtained 
from it by + ^ — ]. What is the metric solving Einstein equations (|2]) at x^,x^ > 0, 
i.e., in the forward light-cone region, as shown in Fig. [Sp 

To date there exists no analytical solution of this problem. On the analytical side 
the problem was tackled perturbatively order-by-order in the EMT's of both shock waves 
[25| |26] and by an all-order resummation of the perturbation theory in the EMT of one 
of the shocks [27] . One usually approximates the nuclear profiles by delta- functions: 
defining ti{x-) = ^ (Ti++(x-)) and t2{x+) = ^ (T2~(a;+)) for the EMT's of the 
two shock waves, we have 

ti{x-) = fli 6{x^), t2(x+) = fl2 S{x+) (19) 

with [261 [28] 

/il~J9+A?Ai/^ /i2~j92 A^A^/l (20) 

Here Pi, P2 are the large light-cone momenta per nucleon in the two nuclei with atomic 
numbers Ai and A2, while Ai and A2 are the typical transverse momentum scales 
describing the two nuclei. 

The metric in the forward light-cone can be represented as a series in the powers 
of yUi and fi2, corresponding to dimensionless Lorentz- invariant expansion parameters 
/ii (x~)^2;+ and fi2 {x'^)'^ x~ . The idea of solving Einstein equations by an expansion 
in powers of /ii and /i2 was formulated in [25] for the delta-function profiles (fT9|) and 
extended to other more realistic longitudinal profiles in [2B] and to non-trivial transverse 
profiles in [29] . The all-order resummation of the powers of /i2 keeping fii at the lowest 
non-trivial order was carried out in [27] resulting in the following expectation value of 
the EMT in the forward light-cone (only non-zero components are listed): 
/T++\ - 4/ii/^2 {x+Y9{x+)9{x-) 

{T--) = ^J{x+)e{x-' ^ 



2vr2 ^ ' ' '2ii2{x+Y [l + 8/X2(x+)2x-]'/' 



X 3-3y^l + 8/i2(a;+)2x- +4/i2(x+)2a;- 
X f 9 + 16 /i2 {x-^f X' -QJl + Sii2 {x+Y X- 



(216) 

- (T-) ^ (T-) ^ ^ 8,,,2X^x-.(x^)^(x-) ^ 

2 7r2 [1 + 8/i2(x+)2a;-] / 



AdS/CFT applications to relativistic heavy ion collisions: a brief review 



9 




Figure 4. Energy density in the shock wave cohision obtained by the numerical 
solution of the problem carried out in |30]. In the Eddington-Finkelstein coordinates 
used in |30| . v is time, while z is the longitudinal coordinate (the collision axis), fi 
is a parameter with dimension of a mass, which is related but is not equal to /ii,/^2 
introduced in the text. (Reprinted figure with permission from [30 . Copyright 2011 
by the American Physical Society.) 



which, by construction, is apphcable in the region defined by fii {x )^ x^ ^ 1. Due to 
this hmitation, the EMT in Eqs. fl2Tal) . fl2T6ll . and fl2Tcjl is not vahd at asymptotically 
late proper times r = when the system is likely to thermalize and should 

be described by ideal hydrodynamics: it does not yield us a new solution of ideal 
hydrodynamics equations. (For instance, < in Eq. (12 lap , which is impossible in 

ideal hydrodynamics.) Rather Eqs. (12 lap . <\21b\i . and (I21cp describe a non-equilibrium 
medium on its way to equilibration. 

The all-order resummation of ^2 {x^)"^ x~ in [27] allowed for one more insight about 
shock wave collisions: smearing the profile of the first shock wave to generate a finite 
width ai 

t^^x-) = ^e{x-)e{ai-x-) (22) 
ai 

one can follow the evolution of the shock wave after the collision, which yields [27] 

(r++) = ^ ^ ^ , for < X- < ai, x+ > 0. (23) 

27r^ ai y'l + 8/i2(x+)2x- 

We see that — t- with increasing x"*". This implies that on the time scales of the 
order of xj^p ~ ^^^^^ the shock wave would lose much of its light cone momentum 
and would "stop", i.e., it would either deviate from its light-cone trajectory or would 
simply dissolve into the produced medium. This "stopping" effect may be due to the 
following mechanism: the partons in the colliding nuclei branch unchecked due to the 
large value of the coupling, which results in the shock waves consisting of a huge number 
of very low-momentum partons. When two systems of low-momentum partons collide, 
the partons easily slow each other down, resulting in the "stopping" of Eq. f l23|) . This 
"stopping" conclusion, unfortunately, is a strong argument against using AdS/CFT to 
describe the whole heavy ion collision, since no such "stopping" behavior appears to 
have been observed experimentally at RHIC and LHC 
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A numerical solution of the shock wave collision problem was carried out in [3U] for 
the shock waves with Gaussian profiles (see also pi]). The energy density obtained in 
[30] is plotted in Fig. HI The results of [30] appear (at least qualitatively) to confirm the 
above analytical conclusion about the fast depletion of the EMT along the shock wave's 
light-cone ( |23l) . In addition, the numerical simulation suggests a rather fast onset of 
viscous hydrodynamics, indicating a fast equilibration of the produced medium. 



3.2. Trapped surface analysis and thermalization 

Interestingly enough, one does not have to solve Einstein equations explicitly to figure 
our the fate of the colliding shock waves system. It is possible to determine whether the 
system reaches thermal equilibrium, that is whether a black hole is created in the bulk, 
by performing a trapped surface analysis [32] [33] . According to the Hawking- Penrose 
theorem, existence of a trapped surface implies that a gravitational collapse (and, hence, 
a black hole creation) is inevitable. A trapped surface analysis for colliding shock waves 
in AdSs was carried out in [311 ESI EEl EZ], with the original work on the subject [31] 
dealing with shock waves with non-trivial transverse coordinate profiles, obtained by 
putting ultrarelativistic point sources in the bulk. We now know that trapped surfaces 
are formed in collisions of a variety of shock waves with and without sources in the bulk 
[M| l35| [36| l37] , indicating that thermalization does take place in heavy ion collisions at 
strong coupling. 

The area of the trapped surface gives one the lower bound on the entropy of the 
black hole, and hence, on the entropy of the produced matter. For the sourceless delta- 
function shock waves of Eqs. ( lT8i) . ( fT9l) the produced entropy per unit transverse area 
A I was found to be 



£^^(2....)-. (24) 

The thermalization (black hole formation) proper time is given by the only dimensionful 
boost-invariant parameter in the problem (/ii fi2) [25| |37] as 



which, as can be seen from Eq. fl20|) . is very short parametrically, being suppressed by 
a power of the center of mass energy s as /xi /i2 ^ Vi ~ ■s. This conclusion was 
supported by the numerical simulations of [30] performed after [251 [37]. Unfortunately, 
substituting realistic values for /ii,/i2 into Eq. (1251) leads to the thermahzation time 
which is too short to agree with the hydrodynamic analyses of the RHIC and LHC data 
[23 [27]. 

It is worth noting here that the numerical simulations of [30] give 0.35 fm/c as the 
time interval between when the Gaussian shocks used in the simulations start to overlap 
significantly to the onset of viscous hydrodynamics in RHIC kinematics. This number 
is somewhat closer to the time scale on the order of 1 — 2 fm/c used in hydrodynamic 
simulations [381 ESI Sffl E] than the time scale of about Tth ~ 0.1 fm/c resulting from 
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Eq. (123]) evaluated for RHIC kinematics. The difference between the two numbers may 
be due in part to the finite longitudinal width of colliding shock waves used in the 
numerical simulations [30] versus the zero effective width of the delta-function shock 
waves (ITQ!) . However, note that to properly compare thermalization time with the 
hydrodynamic simulations [381 l39l HOl HI] one has to start measuring the proper time 
from r = 0, i.e., from the moment of the complete overlap of the two shocks: such 
procedure, when applied to the results of [30], yields the thermalization time of about 
Tth ~ 0.17 fm/c, much closer to Eq. (1251) and further away from the phenomenological 
times. 

Finally, noting again that yUi /i2 ~ ■§ with s the center of mass energy of the collision, 
we get [31] 

4- s'/'. (26) 

Identifying the produced entropy with the number of degrees of freedom, and hence 
with the number of hadrons produced in the collision, we conclude that the AdS/CFT 
prediction is that the hadron multiplicity should grow as dN/drj ~ s^^^ with rj the 
pseudo-rapidity. Unfortunately the experimental data gives a much smaller power of 
energy, dN/drj ~ s^'^^ [42j, possibly signaling again that the early stages of heavy ion 
collisions (in which the majority of the entropy is produced) are not strongly-coupled. 
Note, however, that modifications of AdSs geometry allow for a better agreement with 
the data Il3]. 



4. Evolution of the produced medium: hydrodynamics and the shear 
viscosity bound 

Let us now change the strategy and assume that the actual collision is either weakly- 
coupled or is described by the physics outside the gauge-gravity duality for some other 
reasons. Strong-coupling physics in this scenario becomes dominant at some later 
(though not necessarily very late) proper time r = -\/2 x+ x~ in the collision. We are 
thus confined to the forward light-cone of Fig. 12] without specific information about the 
past of the system (that is, without the shock waves). The question now is whether we 
can say anything specific about the strong-coupling dynamics of such a medium. 

Assuming again that the system is uniform in the transverse direction, we see that 
the dynamics in the forward light cone in 4d depends only on two variables: the proper 
time and the space-time rapidity rjst = (1/2) ln(x4./x_). Furthermore, let us assume 
for simphcity that the matter distribution is T^st-independent [21]: the assumption is 
justified by the weak rapidity-dependence of the particle multiplicity near mid-rapidity 
in the actual heavy ion experiments. The most general metric in AdSs describing rj^t- 
independent medium in the forward light cone in the Feffermann- Graham coordinates 
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IS 



-e"("'") dT^ + 



dVst + e'^"''^ dxl + dz" 



(27) 



with a(r, z), 6(r, z) and c(r, 2;) some arbitrary functions. One can show that in the r]st- 
and x_L-independent case the most general EMT can be written as 



( <r) 








p{t) 








p(r) 




(28) 







P^{r) ) 

at rist = 0, with the energy density e(r), along with the transverse and longitudinal 
pressures p(r) and Psi^r). 

Let us study the late-time asymptotics of this strongly-coupled medium. Following 
one assumes that at late times the energy density of the medium scales as a 



power of proper time, e ~ r with —4 < A < 0. One can then show that in the 

V = z T^l'^ = const, r — > 00 limit the functions a, b and c depend on the scaling variable 

V only |23j. Solving the resulting Einstein equations and requiring that there are no 
singularities in the bulk, Janik and Peschanski |23] found that A = —4/3 and arrived 
at the following metric (cf. f|T6|l ) 

. 2 



ds' 



z^ 




+ dx ,] + dz 



,(29) 



which is the gravity dual of the celebrated Bjorken hydrodynamics [21]. Note that 
now e(r) = 3p(r) = 3p3(r) = {N^ /27c'^) cq/t^^^, in agreement with [21] (cq is an 
arbitrary dimensionful constant). Since e ~ in this conformal plasma, we conclude 
that T ~ r^^/^, again in agreement with [21]: the temperature falls off with time as the 
system cools. The metric ( |29ll has a horizon at 



o ^ 1/4 

-) r^/l 
eo/ 



(30) 



Comparing with Fig. [T] we see that the horizon of the black hole ( l29l) falls deeper into 
the bulk as the time goes on. 

The important conclusion we draw from the result of [23] is that a strongly-coupled 
rjst- and x_L-independent Af = A SYM medium would invariably end up in a Bjorken 
hydrodynamics state. Hence one can derive the late-time asymptotics of this medium 
without an explicit knowledge of the medium's origin. 

It is also interesting to study the approach to the Bjorken hydrodynamics/ Janik- 
Peschanski metric (!29l) . This was done in [441 145] by expanding the coefficients of the 
metric (1271) around the scaling solution (l29ll for late times, 

1 1 

a(r, z) = a{v) + ai{v) ^ + 02(1-) + ... (31) 



r 



^4/3 



I Note that the bulk z-variable in Feffermann- Graham coordinates is, in general, different from z used 
in, say, Eq. for that particular black hole metric the relation is given by Eq. ([T5|) . 
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with similar expansions for 6(r, z) and c(r, z). The coefficients ai, 02, &i, • • • were found 
by solving Einstein equations perturbatively in 1 / r^/^ with the matching near the AdS 
boundary onto viscous hydrodynamics for the EMT of the gauge theory. The latter 
yields the following EMT ^ 

\ 

p(r) + ff 









+ 











p{r) - ^ 2 



(32) 



3 r 



with rj now denoting shear viscosity. In a conformal medium where e ~ the 
shear viscosity scales as rj ~ ~ 1/r in the Bjorken expansion. We hence write 
rj = {N^ /27r'^) Cq^^ r]o/T with 770 a dimensionless constant. Requiring that the resulting 
metric has no singularities in the bulk yields |15] 

^^ 

18 

such that the shear viscosity is [13 H?] 



(33) 



V 



.3/4 
60 



23/2 33/4 ^2 



(34) 



Comparing the EMT of a SYM plasma with e = {N^/2tt^) eo/r^/^ we read off the 
temperature 

21/2 



T(r) 



1/4 



31/4 vrrV3 
which, when used in Eq. flMl) gives 



TT 



(35) 



(36) 



This important relation between the shear viscosity and temperature of the strongly- 
coupled J\f = 4 SYM plasma was derived originally by Kovtun, Policastro, Son, and 
Starinets |1H1 SHI EDI EH |52] using the Kubo formula to relate rj to the absorption cross 
section of a graviton by the AdSSBH. 

Note that combining Eqs. (EE]) and one gets pl[i9| [5nH5T| [52] 

- = i^. (37) 

s Ait 

where s = S/V^d is the entropy density. Based on the consistency of this result for a 
variety of dual geometries it has been conjectured in [511 [52] that the ratio ([37|) is a 
universal lower bound on 77/5 for any theory (including QCD): 

rj , 1 



> 

s All 



(3^ 



This is known as the Kovtun-Son-Starinets (KSS) bound. It was proven for theories 
with gravity duals in [S2]. The possibility of a violation of this bound was discussed 
in [531 EH [55]. The question of the value of rj/ s in QCD has been extensively studied 
in lattice QCD [SH] and by viscous hydrodynamics simulations for heavy ion collisions 
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ESI Sni SI]- The current estimates place rj/s in QCD (extracted from heavy ion 
colhsions data) very close to the KSS bound [H], though one has to remember that 
QCD is not a conformal theory, and ?7/s is T-dependent in QCD unlike TV = 4 SYM 
[56]. 

Returning to the time evolution of the strongly-coupled rapidity-independent 
medium, note that further higher-order corrections have been calculated in 

[171 157]. Reversing the problem, one can ask about the possible early-time asymptotics 
of this medium. It was shown in |58] that for e ~ only even integer A are allowed as 
r — i- 0. The non- negativity of the energy density demands that A = [581 ES]. In such 
a scenario the energy density e starts out as a constant in time and eventually begins 
to fall off like e ~ r~^/^ as dictated by Bjorken hydrodynamics. [Note however that if 
one expands Eqs. ( 121ap . ( ]21&[) . and (12 1 c\j to the lowest order in r (as = re^''"* / \/2) 
one would obtain e ~ [25j, such that A = 2, though the medium produced in a shock 
wave collision is not boost invariant.] For numerical simulations of rapidity-independent 
(boost-invariant) plasmas see IBT]. 

The weakness of the approach presented in this Section is at the early times: one has 
to specify the metric to be used as the initial condition for Einstein equations at some 
small value of r = tq. The EMT of the particles formed in the early stages of a heavy ion 
collision would only specify gj^J at tq in the expansion of the metric near z = 0. Higher 
order metric coefficients g'^^l, i > 6, have to be specified too at tq: in the AdS/CFT 
dictionary they are likely to be dual to some higher-dimensional operators in A/" = 4 
SYM theory. It appears that to insure the matching between the (presumably) weakly- 
coupled early-time dynamics and the later-time dynamics described by AdS/CFT one 
has to specify an infinite tower of expectation values of operators at tq resulting from 
the former to be used to initiate the latter. To date the problem of how to do this in a 
both rigorous and realistic way remains open. 

5. Jet quenching in a strongly-coupled medium 

Many experimental signals of QGP formation and evolution can be calculated in the 
framework of AdS/CFT duality. Above we have discussed the implications of the 
AdS/CFT correspondence for the hydrodynamic description of heavy ion collisions. In 
the spirit of a brief review, we will only mention one more signal of QGP: jet quenching. 
The idea that jets can lose energy in QGP which would lead to their suppression was 
originally proposed by Bjorken in [62], with the first calculations of jet quenching in the 
perturbative QCD framework carried out in [631 E EH 165] . 

In the AdS/CFT framework the question of jet energy loss was first tackled in [8] 
with the goal of determining a strong-coupling value of the jet quenching parameter q 
originally defined in [7] in the perturbative QCD framework. The purely AdS approach 
to jet quenching was first advocated in [661 EZl EH EHl ITQl ITT] and involved calculating 
the drag force on a heavy quark moving through the SYM plasma. At zero temperature 
a heavy quark is dual in AdS to the endpoint of an open Nambu-Goto string terminating 
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on a D7 (flavor) brane wrapped around S^, with the other end of the string stretching 
to the stack of D3 branes at z = oo [72] . For a finite temperature medium, if we work 
in the coordinate system employed in writing down the metric (jlj), we can only follow 
the string down to the AdSSBH horizon. 







1 - 

Lhhmhh|||| D7 brane ^^^MMB 




/ ^ V 




string 






Y 


black hole 



Figure 5. A sketch of the traihng string configuration dual to a heavy quark 
propagating in the strongly-coupled A/" = 4 SYM plasma. 



Let US find the drag force on a heavy quark moving with a constant speed v through 
a supersymmetric plasma (assume that some external force is applied to keep the quark 
velocity constant). To do this, one first has to find the string configuration in the 
AdSSBH background (jlj) dual to the moving quark. This is done by extremizing the 
Nambu-Goto action of an open sting 

Sng = ~ 2^2 / dr da ^J^^detg^b, Qab = 9mn daX^' d^X^ (39) 

with the boundary condition requiring that the string endpoint attached to the D7 
brane moves with velocity v along the x-axis. Here qmn is the AdSSBH metric (jlj) 
and X^ = X^{T,a) are sting coordinates: they specify the mapping from the string 
world-sheet coordinates cr" = (r, a) with a = 0, 1 to spacetime coordinates x'^ . (Note 
that da = d/do"".) The resulting classical string configuration is 

X^\t = t,a = z) = (t, 0, 0, x{t, z), z) (40) 

with [Ml in ig 



X{t, Z) = vt — 



Zh I f Zh + Z\ TT 

arctan 1 — m 



(41) 



z 2 \Zh — zj 2. 
It is illustrated in Fig. |5l The drag force (momentum change) is given by the momentum 
flow down the string from its endpoint attached to the D7 brane, 
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where p is the 3-momentum of the quark and vr^ is the canonical momentum of the 
string defined by 

„ _ dy-det gbc , . 

Using Eqs. fHT]) . ( HOl) . fj43|) in Eq. fj42|) gives us the drag force on a heavy quark in the 
strongly coupled Af = A SYM plasma f66| 1671168] 

dp TT a/a ^9 p , , 

(it 2 M ^ ^ 

where M is the mass of the heavy quark. (The mass appears since p/M = v/y/l — v'^.) 
Eq. (jH]) was derived for a static medium: it was generalized to the case of Bjorken 
hydrodynamics corresponding to the Janik-Peschanski metric fl2^ in [73] with the 
conclusion that Eq. ( 1^^ still applies in this dynamic case if one replaces T — )■ T(r) 
in it, with T(r) given by Eq. (l35ll . 

The AdS/CFT result for the drag force (jH]) allowed for some interesting 
phenomenology of the heavy quark energy loss at RHIC |71]. However, one has to 
be careful applying this large-coupling result to QCD jet physics: due to asymptotic 
freedom, at least some part of the jet coupling to the medium should be perturbative. On 
top of that, the model considered here involves a quark which is being dragged through 
the medium by an external force which prevents it from slowing down, while in real 
life the hard partons simply plow through the medium loosing energy and momentum. 
Modeling of jets in the medium along these more realistic lines was carried out in [75] . 
The conclusion was that, due to strong coupling, the partons are highly likely to branch 
into more partons, distributing the energy democratically between them. In the end 
one obtains an isotropic distribution of partons, very different from the jet cone of 
perturbative QCD and from the jet cones seen in the actual collider experiments (see 
also [76j). This difference between the "AdS jets" and real life jet gives us another 
argument in favor of being cautious in applying AdS results for jet quenching to the 
real-life heavy ion collisions. 

Applications of the AdS/CFT correspondence to jet quenching are by no means 
limited by the calculation we have just presented. Further research addressed issues 
of studying Mach cones generated by the supersonic jets, energy loss of light quarks 
and gluons, momentum broadening of hard partons traversing the medium, quark and 
gluon stopping distance, dissociation and melting of mesons, along with other topics. A 
proper discussion of these important results lies outside of the scope of this brief review. 
We refer the interested reader to [77] for a detailed review of these subjects. 



6. Conclusions 



Our brief review covers several topics on applications of AdS/CFT correspondence to 
heavy ion collisions, with an emphasis on modeling the actual collisions, stemming 
in part from the author's subjective interests and research experience. Indeed, due 
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to the compact format of this review, many important topics in apphcations of 
AdS/CFT correspondence to heavy ion coUisions had to be omitted. Among those 
topics are two-particle correlations (see |19] for a pedagogical presentation of their 
calculation), electromagnetic probes [78], heavy quark potential at finite temperature 
[ZSl Eni EH ESI IH3] , jet quenching-related observables listed above, and efforts to calculate 
the elliptic flow V2 in AdS [29] . 

We hope that we have convinced the reader that practically any heavy ion 
observable can be modeled using the AdS/CFT correspondence, generating interesting 
qualitative and, sometimes, quantitative insight. In comparing the results of AdS/CFT 
calculations to the actual heavy ion data or to lattice QCD simulations one has to be 
careful to remember the many differences between QCD and JV = 4 SYM theory, though 
one hopes that certain universal features are common for both theories. One may also 
hope that future research would put these differences between QCD and M = 4 SYM 
under quantitative theoretical control. 
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